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We propose an Ashkin- Teller like model for elastic response of DNA molecule to external force and 
torque. The base-stacking interaction is described in a simple and uniform way. We obtain the phase 
diagram of dsDNA, and in particular, the transition from B form to the S state induced by stretching 
and twisting. The elastic response of the ssDNA is presented also in a unified formalism. The close 
relation of dsDNA molecule structure with elastic response is shown clearly. The calculated folding 
angle of the dsDNA molecule is 59.2°. 
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I. INTRODUCTION 



The DNA molecule is the basic genetic material. The 
ability of DNA to pack and fold into chromosomes or 
to serve as a template during transcription and replica- 
tion depends on the particular elastic properties of the 
molecule as modified by local interactions. Attribute to 
the rapid development of the technique of directly manip- 
ulating sing le molecule [il, S, S i, [I, H, measurement of 
elastic properties of single DNA molecule under stretch- 
ing or twisting becomes possible recently. An interesting 
new phenomena revealed by experiments is that there 
exists a new state, called S-state which exhibits exotic 
elastic properties distinguished from the B state. The 
dsDNA molecule can be extended to almost as twice as 
its original length before transits to the new state. There 
is a sharp transition between the B-state and the S-state, 
which shows an evidence of high level of cooperativity of 
molecules in the dsDNA. Authors have argued that there 
must be unknowing dsDNA molecular configuration cor- 
responding to the structural transition. 

References fi 0, H i, 113, El, [11, [H, Ei, [H, El, El, 
ISl . Il9i . ,20,. .2h] provided valuable insights. However, dis- 
crepancy (for example, untwisted or helical) still exists. 
Several models were proposed to fit data of the dsDNA 
molecule response to external force. The freely jointed 
chain (FJC) model treats dsDNA simply as a polymer 
joined freely by rigid units and neglect the base-stacking 
interaction between two adjoint base pairs. It fits the 
data at low force limit well for DNA molecule. However, 
the FJC model can not be used to describe the behavior 
of dsDNA molecule which subjected to stretching and 
twisting simultaneously. It is well-known that twisting 
also plays an important role during the structural transi- 
tion process of dsDNA molecule. From a biological point 
of view, torsional stress is indeed popular in the living 
cell and may strongly influence dsDNA functioning. The 



phenomenological model 01, which takes stretching and 
twisting both into account, describes every state with five 
independent parameters and trends no phase transition 
occurred in the course of stretching and twisting. 

In this paper, we present an Ashkin- Teller like model 
to describe the elastic response of dsDNA molecule in a 
unified framework (various states with the same param- 
eter) and try to find the structural transition mechanism 
of dsDNA molecule. 

II. THE ASHKIN-TELLER FORMALISM 

The FJC model describes dsDNA molecule simply 
as joined rigid units. While stretching one end of the 
molecule along the direction of the strand and fixing an- 
other end, one can deduce a force-extension relation [lij 
as 



(^)=eoth('^V-. 



(1) 



where b is effective unit persistent length and z is dis- 
placement of dsDNA molecule along the strand under ex- 
ternal force /. It is well-known that dsDNA is a double- 
stranded linear biopolymer, each strand is a covalently 
linked chain of nucleotides. Nucleotides come in four va- 
rieties and are composed of three distinct parts: a sugar 
called deoxyribose, a phosphate group, and one of four 
units, adenine, guanine, thymine, and cytosine called 
bases. Each of the four bases attached to nucleotides is 
hydrophobic and capable of hydrogen-bonding with an- 
other complementary base. Each arm of the dsDNA is 
formed, in turn, by a covalent bond between the phos- 
phate group of one nucleotide and the hydroxyl group 
of its neighbor nucleotide's sugar. The bond is often re- 
ferred to as the phosphodiester bond and the two arms 
collectively build up the DNA backbone. This classical 
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FIG. 1: FJC model: freely jointed rigid units. 



B-DNA structure is the basic dsDNA conformation found 
inside living cells. Now, it is generally accepted that con- 
tribution from base-stacking is the main reason help to 
form steady two strand helix configuration, and it also 
determines the dsDNA's local transformation. Unfortu- 
nately, the precise picture about base-stacking is still un- 
known. 

The FJC model treats the dsDNA molecule composed 
of two interacting strands simply as the freely joined rigid 
units and neglects the interaction between two adjoint 
bases (base-stacking). It is suspected reasonably that it 
fails to describe high force extension and can not take the 
torque effect into account. Here, we suggest an Ashkin- 
Teller like model (based on the Ashkin- Teller model 22, 
[2^ in statistical mechanics) to describe in a uniform way 
the elastic response of DNA molecule to both external 
force and torque. 

The Hamiltonian of the Ashkin- Teller like model with 
external force (/) and torque (cr) is of the form 



H 



i 

+ {af + ^){s, + n)l (2) 

where s and r are 1-dimension unit orientation vectors 
(takes values of ±1). K{f,a) describes the strength of 
base-stacking interaction between two adjoint base pairs. 
In the case of K{f,a) = 0, the model reduces to two 
decoupled chains, a and 7 are constant parameters. It 
is obvious that the strength of base-stacking interaction 
depends on structure of DNA molecule as well as external 
conditions in which the DNA is studied. In general case 
of the DNA molecule under stretching and twisting, we 
give a simple phenomenological describing of the base- 
stacking interaction strength as 



K{f, cr) = -{O.la - 0.01)6"-'^+'' - 0.6(a/ + 7) - 2 , (3) 

where a denotes the supercoiling degree defined as 

Lk — Lkq 



FIG. 2: The Ashkin- Teller like model consists of two coupled 
chains. A possible transform of configurations corresponds to 
elastic response to external interaction. 



here is the linking number of the two strands whose 
value is the sum of the twist and the writhe, Lkq is the 
linking number of two strains in natural state. 

The average unit bending persistence length can be 
calculated by the formula 



l = -{s + T)+e 



(s + r) = 



(4) 



where b is the unit physical bending persistence length 
(takes the value of 6 = 0.34nm per base pair), and z 
denotes the partition function of DNA molecule 



z= e 

Si,Ti = ±l 



-pEi 



(3 - l/keT 



(5) 



The total relative extension of dsDNA molecule under 
stretching and twisting can be obtained 



Nl 



N 

To 



(6) 
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In Fig. 3, we present a plot of the relative extension of 
dsDNA molecule vs. stretching and twisting. 

The surfacegraphics demonstrates clearly that the 
Ashkin- Teller like model describes well the elastic re- 
sponse of dsDNA molecule under both stretching and 
twisting simultaneously with same parameters for the B 
and S state when the DNA molecule stretched over its 
contour length. 

Fig. 3 exhibits three distinct DNA states (correspond- 
ing to different elastic property) and transitions among 
them induced by external stretching. Furthermore, it re- 
veals quantitatively that the transition between B-form 
to S-state is influenced greatly by external twisting. With 



3 



Relative extension 




It is in good agreement with experimental observations 
0. 

The model gives also the asymmetric behavior of the 
relative extension of the dsDNA molecule between nega- 
tively supercoiled and positively supercoiled region (see 
Fig. 5). At the region of negatively supercoiled, the 
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FIG. 3: Relative extension of dsDNA molecule vs. stretch- 
ing and twisting. The surface-graphics is obtained for a — 
0.04348 and 7 = 0.43478. Other parameters are selected as: 
Lo ~ IS.l^m (from Ref JJ), e = 0.04nm, TV = 71059 and 
b = 0.34nm as physical length per base pair. 




-f=60 
-f=65 

■f=70 

• f=75 
■f=80 
■f=90 

• f=100 



Ttie degree of supercoiling a 



the supercoiling a increasing from negatively to posi- 
tively supercoiled, the transition from the B form to the 
S-state needs larger external force and the cooperative 
transition becomes more sharpened. This is in agreement 
quantitatively with experiments. 
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FIG. 5: Relative extension of the dsDNA molecule under 
twisting with fixed force. The theoretic curves are obtained 
with the same parameters in Fig. 3. 

elastic response of dsDNA behaves differently for differ- 
ent twisting while keeping force a constant. However, 
when cr rises up to about 0.1, all curves get-together and 
have no difference in the course of a increasing continu- 
ously. These results indicate that twisting have much 
more influence on negatively supercoiled dsDNA than 
on positively supercoiled one . The realistic significance 
of the exotic behavior for dsDNA molecule is that, in 
the living, cell dsDNA is usually negatively supercoiled 
(cr K, —0.06). Experiments, 4] really showed asymmet- 
ric behavior of elastic response of dsDNA molecule when 
external force is under 8 pN and the stretched dsDNA 
molecule is smaller than its contour length. Our model 
reveals that the same property still exists after dsDNA 
molecule is stretched over its contour length. We wish 
more elaborate experiments should check the validity of 
the prediction. 



FIG. 4: Relative extension of the dsDNA molecule under 
stretching with fixed twisting. The theoretic curves are ob- 
tained with the same parameters in Fig. 3. 

To further investigate the relative extension of the 
DNA molecule under stretching with fixed twisting and 
to compare with the experiment dataQ, we present the 
force-extension curves with fixed supercoiling a in Fig. 4. 



III. FOLDING ANGLE CALCULATION AND 
REDUCING DESCRIPTION OF SSDNA 

Refs. [H, nil postulated that dsDNA structure should 
be extremely important to its elastic properties and in- 
troduced a structural parameter: the folding angle 
which describes the angle between one of dsDNA back- 
bones and central axis. In the Ashkin- Teller like model, 
value of folding angle can also be obtained 
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FIG. 6: Folding angle vs. relative extension under external 
force(in the case of ct = 0, i.e., torsion free). While the dsDNA 
molecule is stretched to its contour length, the folding angle 
has value tp — 59.2°, which is in good agreement with the 
realistic structure of dsDNA (with ip ~ 62.0°). 



FIG. 7: Relative extension of ssDNA under stretching. The 
theoretic curves are obtained with a = 0.00086 and 7 = 
0.15437. Other parameters are selected as: Lo = 3.9^m (from 
Ref.[^), e = — 0.35nm, A'^ = 3250 and b — 0.6nm as its physi- 
cal segment length. Experiment data are from Ref.Q|. 



(costp) = 
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(7) 



Fig. 6 demonstrates that the elastic behaviors of dsDNA 
molecule under external force are really closely related to 
its structure and the postulation [13, 13| is verified. 

As shown above, the Ashkin- Teller like model describes 
well the elastic behavior of dsDNA molecule. The base- 
stacking interaction between the double stranded struc- 
ture is represented by the strength of "four site" in- 
teraction K{f, a) uniformly. We know that this base- 
stacking interaction is the main factor maintaining the 
stable double-helix structure and resulting in configu- 
rational transformation. Without this interaction the 
model reduced to two decoupled chains. The fact pro- 
vides possibility of describing elastic response of dsDNS 
and ssDNA in a unified frame. 

Fig. 7 really shows a good agreement with ssDNA 
experiment data [^. Other models [3] dealt with the 
elastic response of dsDNA and ssDNA separately and 
have severe deviation in describing of extension for ss- 
DNA when external force is larger than AQOpN. 



IV. CONCLUSIONS 

The proposed Ashkin- Teller like model gives a quan- 
titative description for elastic response of dsDNA under 
both stretching and twisting in a uniform way. It is in 
good agreement with experiment data. We obtained the 
phase diagram of dsDNA, and in particular, the transi- 
tion from B form to the S state induced by stretching and 



twisting. The close relation of dsDNA molecule structure 
with elastic response has been shown clearly. 

Acknowledgments: We thank Yanheng Li and Yue 
Yu for valuable discussions. The work was supported in 
part by NSFC under Grants 10575106 and 10375072. 



* Electronic address: lchangz@m ail.ihep.ac.cn| 
^ Electronic address: pwang@ m ail . ihep . ac . cn | Electronic 
address: yh'zhengQmail.ihep.ac.cn' 
[1] S. B. Smith, L. Finzi, C. Bustamante, Science 258, 1122 
(1992). 

[2] S. B. Smith, Y. Cui, and C. Bustamante, Science 271, 
795 (1996). 

[3] P. Cluzel, A. Lebrun, C. Heller, R. Lavery, J. L. Viovy, 

D. Chatenay, F. Caron, Science 271, 792 (1996). 
[4] T. R. Strick, J. F. AUemand, D. Bensimon, A. Bensimon, 

V. Croquette, Science 271, 1835 (1996). 
[5] M. Rief, H. Clausen- Schaumann, and H. E. Gaub, Nat. 

Struct. Biol. 6, 346 (1999). 
[6] H. Clausen-Schaumann, M. Rief, C. Tolksdorf, and H. E. 

Gaub, Biophys.J. 78,1997 (2000). 
[7] J. F. Leger, G. Romano, A. Sarkar, J. Robert, L. Bour- 

dieu, D. Chatenay, and J. F. Marko, Phys. Rev. Lett. 83, 

1066 (1999). 

[8] A. Lebrun and R. Lavery, Nucl. Acid. Res. 24, 2260 1996. 
[9] J. F. Marko, Europhys.Lett. 38, 183 (1997). 
[10] J. F. Marko and E. D. Siggia, Macromolecules 28, 8759 
(1995). 

[11] R. D. Kamien, T. C. Lubensky, P. Nelson, and C. S. 

O'Hern, Europhys. Lett. 38, 237 (1997). 
[12] C. Bouchiat, M. D. Wang, J. F. AUemand, T. Strick, S. 

M. Block, and V. Croquette, Biophys.J. 76, 409 (1999). 
[13] J. F. Marko, Phys.Rev.E 57, 2134 (1998). 



5 



[14] I. Rouzina and V. A. Bloomfield, Biophys.J. 80, 882 

(2001). 

[15] I. Rouzina and V. A. Bloomfield, Biophys.J. 80, 894 

(2001). 

[16] M. N. Tamashiro and P. Pincus, Phys. Rev. E 63, 021909 

(2001). 

[17] Zhou Haijun, Zhang Yang, and Ou-Yang Zhong-Can, 

Phys.Rcv.Lctt 82, 4560 (1999). 
[18] Zhou Haijun, Zhang Yang, and Ou-Yang Zhong-Can, 

Phys.Rev.E 62, 1045 (2000). 
[19] C. Storm, and P. C. Nelson, Phys.Rev.E 67, 051906-1 



(2003). 

[20] A. Ahsan, J. Rudnick, R. Bminsma Biophys.J. 74,132 
(1998). 

[21] S. Cocco, J. Yan, J. F. Leger, D. Chatenay, and J. F. 

Marko, Phys.Rev.E 70, 011910-1 (2004). 
[22] J. Ashkin and E. Teller, Phys.Rev. 64,178 (1943). 
[23] M. Kohmoto, M. den Nijs, and L. P. Kadanoff, 

Phys.Rev.E 24, 5229 (1981). 
[24] C. Bustamante, J. Marko, E. Siggia, S. Smith, Science. 

265, 1599 (1994). 



